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Abstract 

We have introduced excluded volume effect, which is an important factor to model a realistic 
pedestrian queue, into queueing theory. The probability distributions of pedestrian number 
and pedestrian waiting time in a queue have been calculated exactly. Due to time needed to 
close up the queue, the mean number of pedestrians increases as pedestrian arrival probability 
(A) and leaving probability (p) increase even if the ratio between them (i.e., p = X/fi) remains 
constant. Furthermore, at a given p, the mean waiting time does not increase monotonically 
with the service time (which is inverse to p), a minimum could be reached instead. 

Keywords Queueing Theory, Asymmetric Simple Exclusion Process, Pedestrian Dynamics 
Research Activity Group Applied Integrable Systems 



> 



O 
O 



X 



1. Introduction 

Queueing theory [1] is one of the most famous and im- 
portant theory in these days since it is applied to many 
systems in the real world such as traffic systems [2] , pro- 
duction systems [3] , networks [4] , and so on. The mathe- 
matical formulation for mean waiting time, which is usu- 
ally calculated by Little's theorem [1], is widely used due 
to its simplicity. In the queueing theory, the state of a 
queue is represented by the number of jobs, which are ve- 
hicles, pedestrians, and packets in networks. When there 
are some jobs in the queue, one job is always receiving 
service, and when it leaves the queue, the service for next 
job starts immediately. This phenomenon is suitable for 
a queue of packets since operation for next packet starts 
instantly by a computer. However, it is not realistic for 
a queue of vehicles and pedestrians since there is a delay 
of moving to service window due to the excluded volume 
effect which is not included in the queueing theory. 

The excluded volume effect is studied in detail by ana- 
lyzing the asymmetric simple exclusion process (ASEP) 
[5]. Many traffic models and pedestrian dynamics mod- 
els are developed by extending ASEP [6,7]. They are 
very successful since the excluded volume effect works 
adequately to represent real movement of vehicles and 
pedestrians. Therefore, we introduce the excluded vol- 
ume effect into the queueing theory for the first time in 
this paper to develop a practical theory for a pedestrian 
queue. The probability distributions and the means of 
both pedestrian number and pedestrian waiting time in 
a queue are calculated exactly and compared with those 



obtained from normal queueing models. 

2. Models and Mathematical Analysis 

2. 1 Outline of the Three Queueing Models 

As a comparison, in addition to our excluded vol- 
ume effect introduced queue (E-Queue), we briefly re- 
view normal queue (N-Queue) with continuous time (N- 
Queue (C)) (which is the most famous queueing model 
known as M/M/l [1]), as well as that with discrete time 
(N-Queue (P)). In E-Queue and N-Queue (P), paral- 
lel update is adopted because it is realistic for one di- 
mensional pedestrian dynamics [8]. Fig. 1 is a schematic 
view of N-Queue (P) and E-Queue. At each time step a 
pedestrian arrives at the queue with probability A and a 
pedestrian at the service window (pedestrian A) leaves 
the queue with probability fi. In N-Queue (P) (Fig. 1 
(a)), pedestrian B moves to the service window as soon 
as pedestrian A leaves the queue. In contrast, pedestrian 
B cannot move to the service window at time step t+l in 
E-Queue (Fig. 1 (b)) due to the excluded volume effect, 
he/she moves there at time step t + 2. 

2.2 Master Equations for N-Queue (P) 

Since mathematical analysis on N-Queue (C) is de- 
scribed in detail in many books [1], we start from N- 
Queue (P). The master equations in the stationary state 
are described as follows: 

P(0) = (1 - A)P(0) + (1 - A)/*P(1), (1) 

P(l) = AP(0) + (1 - \)pP{2) 
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Fig. 1. Schematic views of time variation of qucucing states, (a) 
N-Queue (P). (b) E-Queue. The cell at the right end in the 
queue is the service window. A £ [0, 1] and fi £ [0, 1] represent 
the arrival probability and the service probability, respectively. 



+ {A M + (1-A)(1- M )}P(1), (2) 

P(n) = A(l - i£)P{n - 1) + (1 - X)fj,P{n + 1) 

+ {X/j, + (1 — A)(l — n)}P{n) (n>2), (3) 

where P(n) represents the probability that there are n £ 
[0, oo) pedestrians in the queue. Note that the stationary 
state exits only when A < A cr is satisfied. A cr is a critical 
value of A, and when A > A cr , queue length tends to 
infinity. By solving these recurrence equations among 
three terms, we obtain P(n), and P\v{t) {t G [0, oo)) 
(Probability distribution of the waiting time, i.e., time 
between a pedestrian arrives at the queue and leaves 
there.), N — n -P( n ) (Mean number of pedestrians 

in the queue), and W = ^2^L tPw{t) (Mean waiting 
time) are calculated as shown in Tab. 1. 

2.3 Exact Solution for E- Queue 

In E-Queue, the state is determined not only by pedes- 
trian number n. Fortunately, due to deterministic move- 
ment of pedestrians in the queue (i.e., pedestrians move 
one cell in one time step if their proceeding cell is va- 
cant), two consecutive vacant cells never appear in the 
stationary state. As a result, there are 2™ states when 
there are n pedestrians in the queue since we only need 
to consider whether there is a vacant cell or not in front 
of each pedestrian. Schematic views of the stationary 
states in the case n = 2 are depicted in Fig. 2. 

Our target in this paper is to obtain probability dis- 
tributions of pedestrian number and pedestrian waiting 
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Fig. 2. Schematic views of the stationary states of E-Queue in 
the case n = 2. (a) Group A. The service window is occupied by 
a pedestrian, (b) Group B. The service window is vacant. 
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Fig. 3. State transition diagram of E-Queue. 



time, so that the 2" states do not need to be distin- 
guished completely. The important point is that whether 
the service window is occupied or not. Thus, the 2™ 
states are divided into two groups A and B. The ser- 
vice window is occupied in group A and it is vacant in 
group B. For instance, two states belonging to group A, 
and the other two states belonging to group B in the 
case n = 2 as shown in Fig. 2. 

We describe the sum of the probabilities of the station- 
ary states in group A as Pa(^) and that in group B as 
Pb{ti) when there are n pedestrians in the queue. Thus, 
P(n) = PA{n) + P B (n). Note that we have P A (0) = 
and -Pb(O) = P(0). The state transition diagram of E- 
Queue is depicted as Fig. 3 and the master equations in 
the stationary state are described as follows: 

Pa(1) = (1-A)(1- A1 )P A (1) + AP B (0) 

+ (1-A)P S (1), (4) 

P A (n) = X(l-n)P A (n-l) 

+ (1 - - n)P A (n) + \P B (n - 1) 
+ (l-A)Pfl(n) (n>2), (5) 

Pb(0) = (1-A)P B (0) + (1-A) A1 P A (1), (6) 

P B (n) = XnP A {n) + {l-\) f iP A (n + l) 

(»>!)• (7) 

These equations could also be obtained by reducing the 
master equations where all 2" states are distinguished 
[9]. Solving the equations with normalization condition 
P{ n ) = 1, we obtain the solutions in Tab. 1 in the 
case A < A cr . The probability distribution of the waiting 
time P\v(t) is also calculated as 

Q(t,2) 



P w (t)=f(t,l)P(0)+ lf(t-n,n)^P A (n)} + 

n=l 

Q(t-1,2) 

[f(t-n,n + l){{l-fi)P A {n)+P B (n)} 



n=l 



(8) 



where Q(a,b) returns a quotient of a/6, and f(t,n) = 
/U n (l — fiY~ n is the negative binomial distribu- 



n — 1 

tion. The means, i.e., N and W, are also obtained and 
described in Tab. 1. 
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Table 1. Mathematical formulations of physical quantities of three queueing models. The parameter p = A/p represents the ratio 
between the mean service time (1/p) and the mean arrival time (1/A). The expressions of P(n) and Pw(~t) m the table are valid for 
N-Queue (C) when n > and t > 0, respectively, whereas, those are valid in the two parallel update's queues only when n > 1 and 
t > 1. Pw(0) = for N-Queue (P) and E-Queue. Note that the Little's theorem N = XW is satisfied in all three models. 



Type 



N-Queue (C) 



N-Queue (P) 



E-Queue 



Other Expressions for E-Qucuc 



P(0) 



P(n) 



N 



W 



l+M 



1-P 



1 - ■ 



1 — A 



P + Ap 



1-P 



(1-A) 



P w (t) - p) exp[-//(l - p)t] 



M 1 ~ P) f 1 - M 



1-A Vl-A 
1-A p 



1 - A ) Vl-A 



A 1 



(■-A) 



P A (n)=r»- 1 ? — ^— P(0) 
(1 - A)p 



P B (n) = r"- 1 /i A ' P(0) 



(l-A)V 



1 — p + Ap 
(1-A) 2 ' 



3. N-Queue V.S. E-Queue 

In this section, we compare physical quantities of the 
three queueing models. Note that A and p are probabili- 
ties in N-Queue (P) and E-Queue, while they are rates in 
N-Queue (C). Besides, t is discrete in the former models, 
whereas it is continuous in the latter one. However, if we 
regard one time step in the former two models as unit 
time in N-Queue (C), we can compare the physical quan- 
tities of the models in the same time scale. The physical 
quantities are described by three parameters, which are 
A, p, and p as in Tab. 1. However, since p = X/p, there 
are only two independent variables; thus, we use p and p 
as independent ones in the following. Then, A becomes 
a function of p and p described as X(p, p) — pp. 

3.1 Critical Value X cr 

In N-Queues, the critical value A cr = p. In E-Queue, 
X cr = pl (1 + p)- Since p/(l + p) < p , the region where 
stationary state exists in p — p space is smaller in E- 
Queue than in N-Queues. Due to time needed to close 
up vacant cells, which equals to one time step, the length 
of E-Queue diverges easier than that of N-Queues. 

3.2 Mean Pedestrian Number in a Queue 

This subsection focuses on mean pedestrian number 
N in a queue, which are denoted as iVi, N2 and A^3 in 
N-Queue (C), N-Queue (P), and E-Qucue, respectively. 
As can be seen from Tab. 1 and Fig. 4 (a), 

(i) At given p, N±, N 2 and N 3 increase monotonically 
with the increase of p, which is defined as p = \/p. 

(ii) At given p, Ni — N 2 — N 3 when p — > (which 
implies A — > 0). With the increase of p, N\ remains 
unchanged, N 2 decreases and N3 increases. 

In N-Qucue (C), in an infinitesimal time interval At, 
the arrival probability and the leaving probability are 
XAt and pAt, which tends to zero. Thus, changing A 
and p implies a rescaling of time interval, which does 
not change probability distribution of pedestrian num- 
ber. Therefore, Ni is independent of p. In N-Qucue (P), 
P(0) is independent of p; however, since the common 

ratio in P(n) ^ ilp^ P^j decreases with the increase of 
p, P(n) becomes narrower and higher as in Fig. 4 (b). 



Consequently, N 2 decreases as p increases. In E-Queue, 
with the increase of p, P(0) dramatically decreases as 
in Fig. 4 (b) and the common ratio r increases when 
p > 1/(2 + p), so that P(n) becomes wider and lower; 
thus, increases. We would also like to explain this 
phenomenon intuitively. A pedestrian takes "the time of 
closing up" plus "the service time" , denoted as an "ex- 
tended service time" , to go through the service window. 
When p increases the extended service time does not 
sufficiently decrease since the time of closing up remains 




Fig. 4. (a) Mean number of pedestrians N against the ratio be- 
tween service and arrival time p. (b) Probability distributions of 
pedestrian number in a queue. 
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as a constant. At given p, the increase of p implies the 
increase of A, hence, N 3 increases as p increases. 

3.3 Mean Waiting Time in a Queue 

Figure 5 shows the variation of W against 1/p, which 
is a mean service time, in the case p is constant. With 
the increase of 1/p, W increases linearly in N-Qucuc 
(C) and quasi-linearly in N-Queuc (P), which coincides 
with our intuition. W also increases quasi-linearly in E- 
Queue when 1/p is large; however, when 1/p is small, 
it surprisingly achieves minimum W m - ln (= p/(l — ^fp) 2 ) 
at l//i m i n (= p/(l — \fp))i increases as 1/p further de- 
creases, and diverges at l/p CI (= p/(l — p)). Since the 
Little's theorem N — \W [1] is satisfied in all three mod- 
els, W = (1/ p)(N/p). At a given p, with the increase of 
1/p, N does not change in N-Qucuc (C) and increases 
in N-Qucue (P), thus W increases in both cases. In con- 
trast, N decreases with the increase of 1/p in E- Queue, 
hence, the minimum could be reached. 

As we have seen in Fig 5, E-Queue becomes similar to 
N-Qucucs when service time is large and different from 
it when service time is small. Thus, it is useful to know 
quantitatively when we should consider the excluded vol- 
ume effect from the perspective of application. In Fig. 
6, the p-p plane is divided by the curve R = 1.1, where 
R is a ratio between W in N-Qucue (C) and E-Qucue 
described as R(p, p) = W^E-Queue/Wisr-Queue (C) In the 
lower-left region in Fig. 6 R < 1.1, and the difference of 
W is not critically large, so that it may be allowed to 
use N-Queue (C) for simple calculation when both p and 
p are small. In contrast, R > 1.1 and the difference is 
crucial in the upper-right region, therefore, the excluded 
volume effect should be considered when both p and p 
are large. Note that Fig. 6 is an example of the dividing 
curve, and it is possible to depict the other curves by de- 
termining R as a different value. Thus, this diagram is 
helpful to know the error quantitatively and make a de- 
cision whether to use N-Qucuc or E-Qucue for designing 
a queueing system for pedestrians. 

4. Conclusion 

In this paper, we have newly introduced the excluded 
volume effect into the queueing model and obtained the 



probability distributions and the means of the number of 
pedestrians and the waiting time exactly. Such physical 
quantities of the new model are compared with those 
of previous models which do not include the excluded 
volume effect. When the service time is large enough, 
all models become similar; however, when the service 
time becomes small, the time of closing up the queue 
dominates and the waiting time surprisingly increases in 
the excluded volume effect introduced queueing model. 
We also obtain the diagram to examine when the effect 
of the excluded volume effect becomes prominent. 

To construct a more realistic model, the movement 
in the queue should be stochastic since pedestrians do 
not homogenously close up the queue. Furthermore, the 
length of a queue should be also calculated exactly, and 
the validity of the model should be verified by real ex- 
periments or observations in the near future. 

Finally, we would like to mention that this work is 
financially supported by JSPS and JST. 
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